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Abstract. A complex Hilbert space of dimension six supports at least three but not more 
than seven mutually unbiased bases. Two computer-aided analytical methods to tighten these 
bounds are reviewed, based on a discretization of parameter space and on Grobner bases. A third 
algorithmic approach is presented: the non-existence of more than three mutually unbiased bases 
in composite dimensions can be decided by a global optimization method known as semidefinite 
programming. The method is used to confirm that the spectral matrix cannot be part of a 
complete set of seven mutually unbiased bases in dimension six. 



1. Introduction 

Two orthonormal bases B° = {\^),j = 1 . . .d} and B l = = 1 • • • d} of C d ,d G N, are 

mutually unbiased (MU) if the modulus of the inner product of any vectors not in the same 
basis is constant, 

|(^|^)| = -L, for all i,j = l...d. (1) 

In this paper, we will examine three computer-aided analytical methods which could be used 
to prove the conjecture that no more than three MU bases exist in the space C 6 pQ. This open 
problem concerning the geometry of the state space of a low-dimensional quantum system has 
attracted considerable interest in recent years; see [2] for a thorough review of the properties 
and uses of MU bases in complex Hilbert spaces of dimension d. 

First, we review and illustrate an approach which is based on discretising the parameter 
space when searching for vectors that constitute MU bases. This simplification results in only a 
finite — but possibly large — number of states which must be checked. The results obtained will be 
exact once rigorous error bounds have been established [3]. Second, we show that the constraints 
defining an MU constellation can be expressed in terms of coupled polynomial equations. Thus, 
constructing a Grobner basis may show that they have no solution. Finally, we propose a new 
approach by writing the problem in a form suitable to use semidefinite programming. This allows 
one to algorithmically identify the (non-) existence of solutions by means of global optimization 
techniques. 

We will illustrate each of the proposed methods by reproducing known results for MU bases 
in low dimensions [4]. For the sake of simplicity, each of the algorithms will be used to prove the 



non-existence of four MU bases in dimension two. All three techniques have yet to be applied 
successfully to the problem in dimension six: the computational resources required to implement 
them were beyond those available. 

Throughout this paper, we will use the concept of MU constellations [5]: the MU constellation 
{d— 1, A, fx, v}d with < A, fj,, v < d— 1, for example, consists of four sets of orthonormal vectors 
containing d — 1,A, //, and v vectors, respectively, with vectors taken from different sets having 
the required overlap l/Vd (cf. Eqs. ([T])). It is sufficient to list at most (d — 1) vectors in each 
set since a suitable d th vector can be constructed from them. 

The non-existence of any MU constellation smaller than {5, 5, 5, 5}6 = {5 4 }6 proves the non- 
existence of four (or more) MU bases in dimension six. The smallest number of vectors sufficient 
to parameterize a MU constellation {d— 1, A, fx, v}<i is given bys = A + /i + z/ — 1 (the underlying 
equivalence relations between MU constellations can be found in Appendix A of [S]). The 
number s effectively determines the number of variables and equations needed to define a MU 
constellation and, a fortiori, the computational resources required to run any algorithm solving 
the problem. For example, the MU constellation {5,3,3, 3}6 = {5,3 3 }6 contains s = 8 free 
vectors depending on 40 real parameters, whereas the full set of seven MU bases in dimension 
six involves s = 29 vectors depending on 145 real parameters. In order to reduce the problem 
as much as possible, it seems imperative to work with smallest possible constellations such as 
{5,3 3 }6 or {5 2 ,4, 1}6 which have not been found by the extensive numerical searches presented 
in 0. 

2. Discretizing the parameter space 

One of the few rigorous results about MU bases in C 6 is due to Jaming et al. [3]. They have 
shown that a particular class of matrices, the so-called Fourier family [6], cannot be a member 
of a quadruple of MU bases. The Fourier family consists of (6 x 6) matrices F(x\, X2) depending 
on two continuous parameters x\ and x<i\ for x\ = X2 = 0, the general expression reduces to the 
standard Fourier matrix Fj^ = uj^ k /\/6, j, k = ... 5, where ui = exp(27rz/6). We now describe 
the idea of this approach which, as outlined in [Tj, could lead to a no-go theorem in composite 
dimensions such as d = 6. 

Let us assume that there exists a MU constellation C of the form {d — 1, A, fi, u}^, we 
will attempt to find a contradiction. The constellation C can be parameterised by phases 
a = (qi, . . . , a Pd ) T , according to Eq. (8) of [5|. The first step is to approximate this MU 
constellation by another set of vectors with components being restricted to A th roots of unity 
only. Such an approximation is easily achieved by dividing the interval [0, 2ir) into N non- 
overlapping intervals Ij = [(2j — l)n/N, (2j + l)n/N) ( mod 2ir), with j = 1, . . . , N. If the phase 
a,k hes in the interval Ij, we approximate it by the mid-point of Ij, — > a^j = 2irjk/N. The 
mapping a — > a thus sends the vectors of a given MU constellation C to a different collection of 
vectors (denoted by C) which do not necessarily satisfy the MU conditions ([T]). Effectively, the 
Pd-dimensional continuous parameter space a has been replaced by a grid consisting of a finite 
number of N Pd states. The accuracy of the approximation of the MU constellation C by C is 
determined by the value of N, and additional flexibility results from partitioning the range of 
each variable at individually. 

The key step forward made by Jaming et al. [3J has been to establish rigorous bounds on 
the errors of the scalar products introduced by the discretization. These bounds allow one to 
conclude that no MU constellation exists if none can be found for a sufficiently accurate discrete 
approximation. Thus, the search for a contradiction has effectively been reduced to checking 
a finite set of constellations which can be done exhaustively, in principle. If N is small, the 
search is easy to perform; however, the error bounds are not tight for small values of N . Finer 
partitions are necessary which, in turn, come at a computational cost since the number of grid 
points is proportional to N Pd . A hierarchical refinement of the procedure may be used to reduce 



the amount of computational resources j7|. 



2.1. Dimension two 

It is not difficult to illustrate this algorithm by searching for the MU constellation {1 4 }2, i-e. a 
set of four MU bases in the space C 2 — which is known not to exist. If the constellation {1 4 }2 
existed, the four vectors defining it could be written in the form 
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V2 V ei/3 



(2) 



depending on p2 = 2 real parameters a(= a\) and /3(= o^). The conditions resulting from Eqs. 
(P) lead to three constraints, 
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and it is obvious that they have no solution. 

To see the discretization procedure in action, we now show the inconsistency of Eqs. ([3]) by 
checking only a finite number of cases. Assume that there are values of a and f3 which give rise to 
a solution of Eqs. ([3]), and write them as a = ay + Aay and /3 = /3ji + A/Jj/, where ay is the iV th 
root of unity closest to a and differing by Aay from it, etc. The maximum distance between the 
actual value of a and an iV th root is tt/N. Using Taylor's theorem for a differentiable function 

m, 

f(x + h) = f(x) + hf'{x + 9{x)h) , for some 6(x) £ (0, 1) , (4) 
the condition | cos(a/2)| = l/\/2 implies that 
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(5) 



where the triangle inequality 



|«2 1 < \ai + 02) < I <xi I + \a 2 \ , 



(6) 



and the relations |Aaj| < tt/N have been used. The other two constraints lead to inequalities 



-- 1. 



for aj in ([5]). Thus, one needs to check 
. . N. For increasing values of N, the 



which are obtained by substituting either fly or aj 
N 2 inequalities corresponding to the indices j,j' 
numbers otj,/3ji, and aj — /3ji must approach odd multiples of tt/2 which is impossible. Thus, 
for a sufficiently large finite value of N, the inequalities have no solutions in terms of N th roots 
of unity and in some intervals centered around them — thus covering the entire parameter space. 



2.2. Dimension six 

Jaming et al. have used this method to exclude the entire Fourier family, F(xi,X2), from a 
complete set of MU bases [3]. This was achieved by discretising the fundamental region of the 
Fourier family using N = 180 and two other complete bases using N' = 19. The restriction to 
consider sets of the form {I, F(x±, X2), £> 2 , £> 3 } with two orthonormal bases B 2 and £> 3 , reduces 
the computational complexity considerably but it is not dictated by an inherent limitation of the 
method. Interestingly, the non-existence of a finite projective plane was shown by an exhaustive 
search [5], Since the existence of a complete set of MU bases shares some properties with the 
existence of finite projective planes [9], discretization appears a promising avenue. 



3. Using Grobner bases 

The second method we review relies on a technique to algorithmically search for solutions of 
coupled polynomial equations. The approach is based on the construction of a Grobner basis [TO] 
for these equations, a tool developed in commutative algebra. It has been applied successfully 
in [11] to show that no two MU bases of a particular form can be supplemented by a third basis 
plus a single vector which is MU to those of the three bases. Further rigorous results excluding 
many candidate bases from a complete set of MU bases have been reported in [T2]. We first 
formulate the approach for arbitrary dimensions d, then work through an explicit example in 
the space C 2 , and finally discuss the prospects of applying this approach in dimension six. 

Any MU constellation of type {d — 1,A,//, i/}^ will correspond to a point on a s(d — 1)- 
dimensional hypertorus, with s = \ + n + v — 1 [5]. In a parameterization analogous to the one 
given in ([2]), only phase factors of modulus one each depending on a single real parameter 
will occur. Upon expressing each phase factor as e ia i = Xj + iyj, supplemented by the 
conditions x 2 + y| = 1, we are led to a parameterization which requires 2s{d — 1) real variables 

X =(x 1 ,x 2 ,...,y 2s(d _ 1) )GR 2s( - d - 1 \ 

The conditions for vectors, orthonormal in sets with d — l,X,fi and v elements, to be MU 
now turn into N multivariate polynomial equations, 

p,(x) = 0, j = l...N, (7) 

with 

N = l -{s + l)(s - 1) + I(A 2 + /i 2 + v 2 ) + s(d - 1) . (8) 

Any solution of the Eqs. (J7|) gives rise to a MU constellation of type {d — 1, A, fj,, v}d- 
Geometrically speaking, we wish to describe the variety 

V = {x G R 2s ^ : Pl (x) = 0, . . . , Pjv (x) = 0} , (9) 

which is a subset of the space M 25 ^ -1 ). No MU constellation of the form {d — 1, A, fj,, u}^ exists 
if and only if the associated variety is empty, V = 0. The set V has an algebraical description 
in terms of an ideal I = (p±, . . . ,pn) which is generated by the polynomials pi(x) to pjy(x). 
This ideal I consists of all linear combinations a of the polynomials Pj(x) with coefficients rj 
polynomial in the variables x, 

s 

o = ^r i (x)p i (x). (10) 

3=1 

It is important to note that the ideal corresponds to the variety over the algebraic closure of the 
coefficient field, here the complex numbers. 

Having re-cast the problem of identifying constellations of MU vectors in terms of ideals, 
it becomes possible to apply methods from commutative algebraic geometry. The equations 
([7]) have no (real or complex-valued) solutions if the polynomial 1 is contained in the ideal 
(pi, . . . ,Pn}- The construction of a Grobner basis G of the ideal I would allow us to check this 
property since G would be given by the set G = {1}. The converse is not necessarily true: the 
equations defining a MU constellation may have no solutions over the real numbers but the ideal 
is non-empty due to the existence of complex solutions. The equation x 2 + 1 = 0, for example, 
has no real solutions but the ideal / = (x 2 + 1) is not generated by {1}. In a fortunate situation, 
the variety is found to be empty over the complex numbers which would then constitute a proof 
that a complete set of MU bases does not exist in dimension six. 



3.1. Dimension two 

We now confirm the validity of the algorithm based on Grobner bases by proving (again) that no 
four MU bases exist in dimension two. First, we write down the system of coupled polynomial 
equations, the solutions to which would define the MU constellation {1 4 }2- Four real variables 
{x%, X2, 2/1, 2/2} parameterize the candidate vectors shown in Eq. d2|), 



1 



1 



1 



supplemented by the constraints 



(11) 



pi(x) = xf + yf-l 
P2(x) 
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(12) 



The constraints defining a MU constellation {1 4 }2 read explicitly 



p 3 (x) = (1 + Xi) 2 + y\ - 2 
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(13) 



1 + xxx 2 + 2/12/2) 2 + (xiy2 - x 2 y\) 



0. 



It is straightforward to derive that these equations have no solution, a fact which, for the purpose 
of illustration, we will confirm by showing that they have G = {1} as a Grobner basis. 

We have constructed the Grobner basis for Eqs. ()12ll3p using a package called FGb [T3] 
implemented in Maple [13] . A standard desktop PC outputs the desired Grobner basis in 0.016 
seconds, and it is indeed equal to all polynomials over the complex numbers, 



(p 1 (x),...,p 6 (x)) = (l). 



(14) 



A direct confirmation of this result follows from writing down the coefficient polynomials used 
to construct the Grobner basis, namely 
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(15) 
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It is not difficult to verify that the following relation holds, 

5 



J^rj(x)pj(x) 



1 



(16) 



which establishes that G = {1}. It follows that the variety V = {x £ M 4 : pi(x) = 0, . . . ,ps(x) = 
0} is empty implying that no four MU bases exist in the space C 2 . 



3.2. Dimension six 

Unfortunately, 16GB of memory have not been sufficient to decide if the Grobner bases 
corresponding to the equations generated by the constellations {5,3 3 }6 or {5 2 ,4, 1}6 contain 
the element 1. The computations run out of memory before the algorithm terminates 
since the number of variables and equations is considerable, even for these smallest possible 
interesting constellations in dimensions six. For example, the constellation {5 2 ,4, 1}6 requires 
the construction of a Grobner basis of an ideal generated by 61 equations of degree 4 in 90 real 
variables. 

4. Using semidefinite programming 

Semidefinite programming [15] is a powerful tool to obtain rigorous results aided by a computer 
which has been used in areas such as control theory and combinatorial optimization [16] . Recent 
years have seen first applications in quantum information theory to solve problems such as 
deciding whether a given mixed state p is entangled or not [T7] - [20] . The success of semidefinite 
programming stems from the fact that the type of problems studied are efficiently solvable by 
a computer. In addition, there is a duality theorem which provides a certificate allowing one to 
directly check the result. For example, when applied to the separability problem, a semidefinite 
program decides whether a given mixed state is entangled — if it is, it outputs an entanglement 
witness [21], that is, a hyperplane which separates the entangled state at hand from the set of 
all separable states. 

Another interesting application of semidefinite programming to problems about finite 
dimensional Hilbert spaces is to the compatibility problem for tripartite quantum systems. Here 
one asks if there exists a single state of the entire system given the states of all (proper) reduced 
states. Hall [22] has cast the compatibility problem in the form a semidefinite program and used 
it to disprove a conjecture of Butterley et al. [23]. In this case, the solution to the dual problem 
resulted in a certificate, called an incompatibility witness, which proves that a particular set of 
reduced states are not compatible with any multipartite state. 

A semidefinite program (SDP) is an algorithm for solving an optimization problem of the 
form 



where c is a fixed vector, and the decision variables x are constrained by the requirement that 
the matrix F{x) = x\F\ + • • • + x n F n — Fq be positive semidefinite, with symmetric n x n 
matrices Fq,F\, ... ,F n [15]. This is a convex optimization problem. A linear program, where 
the constrains have the form -F(x) = diag(Ax — b), is an example of a SDP with important 
applications in economics [24j or finding optimal network flows |25| . 

A SDP can also be used to solve non-linear problems as long as they are convex. For example, 
the problem 



minimise c x 




(17) 




minimise — — 
a 1 x 

subject to Ax > b 



(18) 



can be re-written as a SDP [15 1 



minimise t 
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diag(Ax - b) 
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(19) 



It is possible to cast the existence problem of MU constellations as a semidefinite program 
|26j . The equations f>j(x) = defining MU constellations are, unfortunately, not convex since 
they involve fourth-order polynomials. However, all is not lost: it is, possible to extend the 
applicability of this method by using tools from non-convex optimization [27, [28]. Lasserre has 
shown that, upon relaxing the non-convex constraints, one can define a hierachy of semidefinite 
programs representing ever better approximations to the original one [29j, at the cost of 
introducing additional decision variables. At each level, this method of relaxations either 
rules out the existence of a MU constellation or is inconclusive requiring an additional step 
of relaxation. Each iteration inevitably leads to a computationally more difficult problem but 
the remarkable work of Lasserre [29] ensures that after a finite number of iterations the exact 
solution will have been obtained. In other words, the hierarchy is asymptotically complete. It is 
this type of approach which has been applied successfully to problems from quantum information 

Let us now show how to express the existence problem of a MU constellation as an 
optimization problem. To do so, 

pick any one of the polynomials from ^7ty, say pi(x), and find the minimum value 
o/p|(x) subject to the condition that the variables satisfy the remaining constraints, 
p 2 (x) = 0, . . . ,pjv(x) = 0. 

Effectively, we seek to minimise the value of pf(x) in the solution space of the remaining 
polynomials. Lasserre's method of relaxations then allows one to find a lower bound Bl(t) 
for the function p\ (x), where r is the degree of the relaxation. If for some degree of relaxation, 
r = 2,3, .. . one finds a positive bound, Bl(t) > 0, then no MU constellation of the desired 
type exists. The dual formulation of the SDP would automatically provide a certificate of non- 
existence, offering an independent way of verifying the result. However, the original problem 
is transformed by the relaxation steps and so this "witness of non-existence" is less intuitively 
connected to the existence of an MU constellation than, for example, an entanglement witness. 

Here is a sketch of the algorithm which is capable to determine whether a MU constellation 
{d — 1, A, n.v}d exists: 

(i) write down the polynomial equations which define the MU constellation; 

(ii) generate a SDP at the lowest possible level of relaxation (r = 2); 

(iii) solve the resulting SDP — if a positive global lower bound Bl(t) is found then no MU 
constellation {d — l,A,/i.z/}d exists; otherwise, repeat Steps (ii) and (iii) at the next level 
of relaxation, r := r + 1. 

If no MU constellation exists, the algorithm is guaranteed to find a positive lower bound. As r 
increases, the global lower bounds B^{r) converge monotonically to the exact global minimum 
of the function, B°^ 1 . If the sought-after MU constellation does exist, the algorithm will find an 
explicit parameterisation thereof, with a high level of numerical accuracy. 

4-1- Dimension two 

To begin, we now spell out and run a SDP algorithm to show once more that that there are no 
more than three MU bases in dimension two. More explicitly, we will calculate a strictly positive 
global lower bound Bi(r) for the polynomial 



p?(x) = (x? + y f-l) 2 , 



(20) 



subject to the requirement that the remaining polynomials in Eqs. (|12p and (|13|) vanish (Step 
(i)). This leaves us with the following minimization problem, 

min (x\ + y\ — l) 2 
subject to x\ + y\ — 1 = 

(l + xi) 2 + y 2 -2 = (21) 

(1 + x 2 ) 2 + y\ - 2 = 

(1 + xix 2 + yiy2) 2 + (^12/2 - ^2yi) 2 -2 = 

with four decision variables x = (x±, . . . , 2/2)- The MU constellation {1 4 }2 exists if and only if no 
global lower bound of pf(x) exceeds the value zero, or Bi{r) < for all r. Any positive bound 
implies the non-existence of a MU constellation {1 4 }2- 

Using the Matlab package gloptipoly3 [31] which is based on the theory presented in [32] . 
we are able to convert the problem (|21[) into a semidefinite program, thus completing Step (ii). 
The resulting SDP has been solved using the SeDuMi MatLab package developed by Strum et al. 
[33], corresponding to Step (hi). The results of the computations at three levels of relaxation, 
r = 2,3,4, are presented in Table [TJ We find that already the first level of relaxation, r = 2, 
provides a positive lower bound for the polynomial p\ (x), B^{2) > 0. This result confirms that 
semidefinite programming is capable to positively identify non-existing MU constellations. 

The third and fourth columns of Table [1] show how the size of the SDP grows as we increase 
the level of relaxation. The number of decision variables, Nd, jumps from 69 to 494. Similarly, 
the dimension of matrices defining the semidefinite constraint -F(x) > is almost five times 
larger at r = 4 than at the lowest level of relaxation, r = 2. The algorithm proves to be very 
efficient taking only 0.11 seconds on a desktop PC to convert the original problem and solve the 
SDP for r = 2. This time rises to 1.67 seconds when r = 4. 
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Table 1. Lower bounds of the minimization problem defined in Eq. (|21|) . with levels of 
relaxation denoted by r, while and F give the number of decision variables and the size 
of the semidefinite inequalities in the resulting SDP, respectively. 

The high level of numerical accuracy achieved by the optimization program SeDuMi allows us, 
in fact, to find an analytic expression for the lower bound, namely B°^ 1 = (1 — \/3) 2 - At the third 
and fourth levels of relaxation, the lower bound is optimal in that the function pf (x) reaches the 
value Bl. It is possible to output the parameter values x which achieve this lower bound. There 
are two sets of vectors corresponding to the global minimum value B° L pt = (1 - ^3) 2 given by 



where a = (v3 — l)/2- Interestingly, the sets V+ and V- both contain three MU bases plus one 
additional vector. 



4-2. Dimension six 

According to Eq. ([8]) the number of variables parameterising candidates for MU constellation 
increases quadratically with the number of undetermined vectors. Thus, a successful application 
of semidefinite programming in dimension two does not guarantee that a similar algorithm will 
terminate in dimension six. As for Grobner bases, finite computational resources may limit the 
type of problems which can be solved. We will consider two MU constellations in C 6 with small 
and large numbers of variables, respectively. 
The columns of the spectral matrix [33] 
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define an orthonormal basis of C 6 since S is unitary. It has been shown that the pair {/, S} 
cannot be extended to a triple of MU bases [12] (the identity / is associated with the standard 
basis of C 6 ). Furthermore, there is no pair of orthogonal vectors {\u), \v)} MU to the column 
vectors of I and S. The two undetermined vectors {\u), \v)} depend on 20 real variables which 
must satisfy 21 constraints given by 4 th order polynomials. Due to the relatively small number 
of variables involved, one may expect that this property can be checked using semidefinite 
programming. Upon transforming the original problem, the equivalent SDP has been found 
to require 10, 625 decision variables satisfying 4, 851 linear constraints while the semidefinite 
inequalities involve matrices of dimension 231 x 231 — at the lowest level of relaxation, r = 2. It 
took approximately 210 minutes and 5.4G of memory for the SDP to terminate, resulting in a 
global lower bound for the square of one of the polynomials defined by the MU conditions, 

B L {2) = 2.28 x 10~ 8 . (24) 

Since -Bz,(2) is positive, the SDP algorithm confirms the result of [12] : no two orthonormal 
vectors can be found which are mutually unbiased to the column vectors of the pair {/, S}. At 
the next level of relaxation, r = 3, the resulting SDP is already too large to terminate: there are 
230,229 variables while the constraints involve semidefinite matrices of dimension 1771 x 1771. 
Consequently, we have not been able to improve the lower bound (I24p . 

It should not come as a surprise our attempts to construct the SDP for the constellation 
{5 3 , 1}6 have not been successful, even at the lowest level of relaxation. Judging by the increase 
in the number of decision variables and the size of the semidefinite constraints seen in the 
previous examples, it is likely to be very large. Solving the resulting SDP would seem optimistic. 
However, it is possible that the resulting SDP has some specific structure which would allow 
one to use appropriate techniques. 

5. Summary and outlook 

We have reviewed two known computer-aided methods suitable to obtain rigorous bounds on 
the number of MU bases existing in the space C d , based on a discretization of parameter space 
and on Grobner bases. Both methods are currently too costly from a computational point of 
view to decide on the (non-) existence of more than three MU bases in dimension six. 

In addition, we have rephrased the problem in terms of semidefinite programming which 
allows one to use rigorous methods of global optimization. This approach has been shown to 
reproduce two known results, the non-existence of four MU bases in dimension two, and the fact 



that the spectral matrix S cannot be part of a complete set of seven mutually unbiased bases 
in dimension six. 

Seemingly innocent questions about sets of mutually unbiased vectors in low-dimensional 
Hilbert spaces such as C 6 continue to resist a variety of rigorous attempts of solution, confirming 
once more the dictum [35] 

Hilbert space is a big place. 

Interestingly some MU constellations in dimension six are just within reach of numerical 
methods [5] which, however, only make their non-existence plausible. It remains to be seen 
whether a proof of the (non-) existence of a MU constellation such as {5 3 , 1}6 will be found 
first by a computer-aided approach or by a purely analytic reasoning. It is safe to say though 
that obtaining a proof will be a matter of time only: the algorithmic methods presented here 
ascertain that we do not deal with an undecidable question when searching for complete sets of 
MU bases in composite dimensions. 
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